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ABSTRACT 

We analyzed simultaneous archival XMM-Newton and RXTE observations of the X- 
ray binary and black hole candidate Swift J1753.5— 0127. In a previous analysis of 
the same data a soft thermal component was found in the X-ray spectrum, and the 
presence of an accretion disk extending close to the innermost stable circular orbit was 
proposed. This is in contrast with the standard picture in which the accretion disk is 
truncated at large radii in the low/hard state. We tested a number of spectral models 
and we found that several of them fit the observed spectra without the need of a soft 
disk- like component. This result implies that the classical paradigm of a truncated 
accretion disk in the low/hard state can not be ruled out by these data. We further 
discovered a broad iron emission line between 6 and 7 keV in these data. From fits to 
the line profile we found an inner disk radius that ranges between ~6-16 gravitational 
radii, which can be in fact much larger, up to ~250 gravitational radii, depending on 
the model used to fit the continuum and the line. We discuss the implications of these 
results in the context of a fully or partially truncated accretion disk. 

Key words: accretion disks, disk states - stars: individual (Swift J1753.5— 0127) - 
X-rays: binaries 



1 INTRODUCTION 

Swift J1753. 5-0127 (hereafter Swift J1753) was discovered 
with the Swift Burst Alert Telescope on June 30 2005 as a 
hard X-ray source (E > 15 keV; Palmer et al. 2005). Swift 
J1753 was also detected at energies below 10 keV with the 
Swift X-Ray Telescope (XRT), and in the UV with the Swift 
UV/Optical Telescope (Morris et al. 2005; Still et al. 2005). 
Using optical spectroscopy, Torres et al. (2005) found double 
peaked Ha emission lines, which indicate that the system is 
a low-mass X-ray binary. The source was also detected in 
radio, which might evidence jet activity (Fender et al. 2005; 
Cadolle Bel et al. 2006). Based on the hard spectral behav- 
ior, Cadolle Bel et al. (2005) proposed that Swift J1753 is a 
black hole candidate. Although, it is generally assumed that 
Swift J1753 contains a black hole, there is yet no dynamical 
confirmation of its nature. Bearing this in mind, in the rest 
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of this paper we treat Swift J 1753 as a black hole candidate 
(BHC). 

Rossi X-ray Timing Explorer (RXTE) and INTE- 
GRAL observations showed that Swift J 1753 never left the 
low/hard state (LHS) during the whole outburst (Cadolle 
Bel et al. 2007; Zhang et al. 2007). The LHS of X-ray bi- 
naries is characterized by a hard power-law X-ray spectrum 
(photon index, F~1.7; e.g., Mendez & van der Klis 1997; 
Remillard & McClintock 2006) , which is usually interpreted 
as the result of Comptonization of thermal photons from an 
optically thick accretion disk by hot electrons (e.g., Esin et 
al. 1997). Observations of sources with low absorption pro- 
vide evidence that in the LHS the accretion disk is cool and 
is truncated at a large radius (e.g., McClintock et al. 2001). 
In the last years, a contribution of a jet to the X-ray emis- 
sion in the LHS has been considered (Markoff et al. 2001; 
Fender, Belloni & Gallo 2004). 

During an outburst, the mass accretion rate onto the 
black hole increases and a transition to the high/soft state 
(HSS) can occur (for reviews, see Homan & Belloni 2005; 
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McClintock & Remillard 2006). The high/soft state is dom- 
inated by thermal emission below ~5 keV, likely produced 
by an optically thick and geometrically thin accretion disk 
with an inner radius at, or close to, the innermost stable 
circular orbit (ISCO; Shakura & Sunyaev 1973). Due to the 
increase of thermal photons, the Comptonizing plasma in 
the system cools down and the spectrum of the source be- 
comes soft. A part of the Comptonized photons is reflected 
off the accretion disk, and produces an iron Ka line between 
6 and 7 keV (e.g., Fabian et al. 1989; Laor 1991), where the 
width of the emission line is related to the size of the inner 
disk radius. In this interpretation, if the disk is truncated 
at large radii in the LHS, no or narrow iron-emission lines 
are expected, whereas the iron line would be broad in the 
HSS where the inner disk radius is close to the ISCO and 
the relativistic effects are strong. 

Recent high-resolution spectral data of the BHC X- 
ray binaries Swift J1753, GX 339-4 and XTE J1817-330, 
showed a soft thermal component in the spectra while the 
sources were in the LHS (Miller et al. 2006a, hereafter M06; 
Miller et al. 2006b; Rykoff et al. 2007). It was suggested 
that the soft component is caused by thermal emission from 
an accretion disk which, even in the LHS, extends down to 
the ISCO. The same data set of XTE J1817-330 was re- 
analyzed by Gierlihski et al. (2008), who showed that the 
spectrum could also be fit without the need of a thermal 
component. 

Here we analyzed a simultaneous archival XMM- 
Newton/ RXTE (the same data as reported in M06) ob- 
servation of Swift J1753, and we tested whether a thermal 
disk-like component is required to fit the X-ray spectra. We 
explored a range of continuum models, beyond the models 
described in M06, with and without a disk component. We 
also report on the presence of a broad Fe line in the spectra. 
In Sj2]we describe the extraction and reduction of the data, 
in we present the results of fitting the data with different 
continuum models, and in Sj4] we discuss the results of the 
continuum models and the line components with respect to 
the issue of a radially truncated accretion disk. 



2 DATA REDUCTION 

We analyzed archival observations of Swift J1753 simulta- 
neously taken with XMM-Newton and RXTE on March 24, 
2006, i.e. ~270 days after the start of the outburst. XMM- 
Newton started its ~42 ks observation at 16:00:31 (UT), and 
the RXTE observation started at 17:32:16 (UT). In Fig. [T] 
we show the date of this observation on top of the All Sky 
Monitor (ASM) RXTE light curve of Swift J1753, starting 
from its outburst in 2005. In the next subsections we de- 
scribe the procedure that we followed to extract the spectra 
for each of the different telescopes. 

2.1 XMM-Newton 

We reduced the XMM-Newton Observation Data Files 
(ODF) using SAS version 7.1.2 and we applied the latest 
calibration files. We extracted the event files for the EPIC- 
pn (hereafter PN) and both EPIC-MOS cameras (MOS1 
and MOS2) using the tasks epproc and emproc, respectively. 
We applied standard filtering and examined the light curves 



for background flares. No flares are present and we used 
the whole exposure for our analysis. We checked the filtered 
event files for photon pile-up by running the task epatplot. 
No pile-up is apparent in the PN and MOS1 data since both 
cameras were operated in 'timing' mode. The MOS1 data, 
however, show a sharp drop in the histogram in which counts 
versus RAWX3 are plotted. This feature is due to a physi- 
cal damage of some pixels in the MOS1 camera as a conse- 
quence of a micro-meteorite impact. For this reason, we ex- 
cluded the MOS1 data from our analysis. The MOS2 camera 
was operated in 'full-frame' mode and suffered from mod- 
erate photon pile-up. We therefore excluded the central 22" 
region of the source point spread function (PSF) and con- 
firmed, using the task epatplot, that the remaining MOS2 
data are not affected by pile-up. The ODF files also include 
high resolution spectra at energies below 2 keV, taken with 
the Reflection Grating Spectrometer (RGS). Since the PN 
and MOS cameras cover most of the range of the RGS spec- 
tra, we do not use the RGS data in this paper. 

We extracted source and background spectra applying 
the standard procedures. The source spectrum for the PN 
data is extracted using RAWX in [28:48], and the back- 
ground in [5:25]. For MOS2, we generated the source spec- 
trum using an annulus with inner radius of 22" and outer 
radius of 68" centered on the source location. We extracted 
the MOS2 background spectrum from a circular region with 
a radius of 160" located 10.5' away from the position of the 
source. We created photon redistribution matrices and an- 
cillary files for the source spectra using the SAS tools rmf gen 
and arf gen, respectively. 

We rebinned the source spectra using the tool pharbr0 
(M. Guainazzi, private communication), such that the num- 
ber of bins per resolution element of the PN and MOS2 
spectra is 3. We checked that after rebinning, all channels 
have at least 15 counts in the energy range used in our fits. 
We note that differently from us, M06 grouped the PN data 
to have at least 10 counts per bin, and they therefore over- 
sampled the intrinsic resolution of the PN data by a factor 
of ~8. The cross calibration between PN and MOS agrees 
to 7% from 0.4 keV to 12 keV0, and we found that only be- 
low 2 keV the systematic difference in the relative effective 
areas of these cameras went up to 7%. We therefore added 
a 7% systematical error to the PN and MOS2 spectra below 
2 keV. 

2.2 RXTE 

We reduced the RXTE Proportional Counter Array (PCA) 
and High Energy X-ray Timing Experiment (HEXTE) data 
using the HEASOFT tools version 6.4. We applied stan- 
dard filtering and obtained PCA and HEXTE exposures 

1 In timing mode, one of the spatial coordinates of the CCD is 
compressed in order to increase the time resolution. For XMM- 
Newton the spatial coordinate that is not compressed in this mode 
is called RAWX. 

2 The pharbn tool allows to rebin the energy channels of the spec- 
tra for different instruments taking into account the instrument's 
spectral resolution. It allows to rebin both on minimum number 
of bins per resolution element and on minimum number of counts 
per bin 

3 http:/ /xmm2.esac.esa.int/docs/documents/CAL- TN-0018.pdf 
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Figure 1. RXTE /ASM light curve of Swift J1753 with daily points, starting at the outburst in 2005 (MJD 53551). The XMM- 
Newton/RXTE observation, ~270 days after the start of the outburst, is represented by the vertical bar. Note that 3 years after the 
start of the outburst, the source is still active. 



which are ~2.3 ks and ~0.8 ks, respectively. Using the tool 
saextrct, we extracted PC A spectra from PCU2 only, be- 
ing this the best calibrated detector and the only one which 
is always on. For the HEXTE data we generated the spectra 
using cluster-B events only, since after January 2006 cluster 
A stopped rocking and could no longer measure the back- 
ground. We estimated the PCA background and measured 
the HEXTE background using the standard RXTE tools 
pcabackest and hxtback, respectively. We further built in- 
strument response files for the PCA and HEXTE data using 
pcarsp and hxtrsp, respectively. The energy channels of the 
RXTE spectra have more than 20 counts per energy bin so 
that x 2 statistics can be applied, therefore no channel rebin- 
ning is required. Following M06, we added a 0.6% systematic 
error to the PCA spectra using the FTOOL grppha, and we 
did not add any systematic error to the HEXTE data. 



3 DATA ANALYSIS AND RESULTS 

For our spectral analysis of Swift J1753 we used XSPEC 
vll (Arnaud 1996). We simultaneously fit the spectra of the 
different instruments, where we restricted the PN and the 
MOS2 spectra to an energy range of 0.6-10 keV and 0.6-9.4 
keV, respectively. For the PCA and HEXTE spectra we used 
the energy ranges 3-20 keV and 20-100 keV, respectively. 
To investigate the nature of the putative disk emission in 
the LHS (as reported in M06), we tested several models to 
describe the continuum spectrum of Swift J1753. 

3.1 General fit procedure 

Each of the spectral continuum models we present below in- 
cludes the effect of interstellar absorption in the direction of 
Swift J1753, using the phabs component with cross-sections 
from Baluciriska- Church & McCammon (1992) and abun- 
dances from Anders & Grevesse (1989). The total column 
density in the direction of the source, as derived from HI 
data, is Nn = 1.7 X 10 21 cm" 2 (Dickey & Lockman, 1990). 
Since this column density provides an upper limit to the in- 
terstellar absorption to Swift J1753, we left A^h free to vary 
during our fits. For an overall normalization between the dif- 
ferent instruments, we multiplied the continuum model by a 
CONSTANT component. All model parameters, except these 



multiplicative constants, were linked between the different 
instruments. 

In the fit residuals of the PN and MOS2 data, two nar- 
row absorption lines near 1.8 keV and 2.2 keV were present. 
This is most likely due to a mismatch in the calibration of 
the edges in the EPIC instruments. To reduce the impact 
of these features in our fit results, we fit Gaussian absorp- 
tion lines, which we found to be at an energy of 1.82±0.01 
keV and 2.25±0.03 keV, with widths of 0.23±0.15 eV and 
2. Oil. 9 eV, respectively. 

There appear to be a broad emission line between 6 and 
7 keV in the XMM-Newton/RXTE spectra. (This emission 
line is best seen in the PN spectra; MOS2 effective area is 
a factor of ~2 smaller than that of PN, and in this case 
the effective area has been reduced due to the excision of 
the central part of the source PSF to avoid pile-up. On the 
other hand, the spectral resolution of RXTE is ~1 keV at 
6 keV.) Similar emission lines, probably due to iron, have 
been reported in other BHC (e.g., Miller et al. 2002b) and 
active galactic nuclei (e.g., Tanaka et al. 1995). We initially 
added a Gaussian component with a variable line width to 
the spectrum to fit this line, but there always remained resid- 
uals around 7 keV which we could only eliminate by adding 
at least another Gaussian to our model. On the contrary, 
a Laor profile (Laor 1991) yields a good fit of the line and 
therefore, in the rest of the paper we use this relativistic line 
profile whenever we fit a line to the spectra, although we also 
tried different relativistic line profiles as described in H3.8I 
For the Laor component we allowed all the parameters to 
be free to vary, except the outer disk radius, which we fixed 
to the default value of 400 R g (R g is the gravitational ra- 
dius of the black hole, defined as R g = GM/c 2 , with G and 
c the common physical constants, and M the mass of the 
black hole). Further, we constrained the line energy to the 
range between 6.4 keV and 7 keV, which is the energy range 
spanned by the Ka lines from neutral to hydrogen-like iron. 

For all fits of the different continuum models, the values 
of the parameters and the 90% confidence errors are given 
in Tab. [1] In the text we only give approximate numbers of 
the parameters. 
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3.2 A power-law model 

The first continuum model we fit was a power-law model 
(pl). This model yields a poor fit (x 2 /i> = 845/502). Fol- 
lowing M06, we added a disk blackbody component diskbb 
(Mitsuda et al. 1984) with a best-fit inner disk tempera- 
ture of ~0.4 keV; the reduced \ 2 m this case is 1.21 (x 2 jv 
are given in Tab. [TJ. The spectra and fit residuals of the 
pl+diskbb model are shown in the left panel of Fig. [2] A 
broad emission line between 6-7 keV is apparent in this 
figure. We added a Laor line profile to the pl+diskbb 
model and the fit improved significantly, with a reduced x 2 
of 1.14. To check whether a fit to the spectrum with only 
a line and a power-law component provides a good fit, we 
fit a pl+Laor model, without the diskbb component. The 
pl+Laor model gives a reduced \ 2 °f 1-30. In the best-fit 
model (i.e. the pl+Laor+diskbb model), the line is re- 
quired at high confidence as indicated by the fact that the 
normalization of the line is 10cr different from zero. The 
line in the best-fit model has an equivalent width of 90.3 
eV, and the diskbb component is required at a confidence 
level of more than 8a. The flux of the diskbb component 
is 7.2xl0~ 12 erg cm" 2 s" 1 (0.6-10 keV), where the unab- 
sorbed flux of the best-fit model is 3.9xl0 -10 erg cm -2 s _1 
(0.6-10 keV). We note that the diskbb normalization that 
we obtained is much smaller than the value that M06 found, 
and the contribution from the disk to the total flux that we 
found is a factor of ~3 less than reported in M06. 

In the case when the distance, inclination angle of the 
disk, and mass of the black hole are known, the inner disk 
radius, J?i n can be obtained from the diskbb normalization 
(No)- However, for the case of Swift J1753 these parameters 
are currently not known. Assuming that the source is closer 
than 10 kpc and has an inclination lower than 63°, since 
eclipses are not yet reported, 7?i n < 1.48 </Nu km. Assuming 
a black hole mass >3 Mq, R ln <0.33 y/~Nn R g . From the 
diskbb normalization we found a non-truncated disk, while 
inferred from the line profile we found Ri n ~15 R g . For an 
inclination angle > 85° , as determined from the line profile, 
Rin can become much larger. 



3.3 Cut-off power-law model 

Next, we added a high-energy cut off to the power-law com- 
ponent. In the case when we let all the parameters free, the 
high-energy cut off became larger than 200 keV. Since this 
value is outside the energy range of the spectra, this com- 
ponent is effectively the same as the simple power law. We 
therefore did not explore this model further. 



3.4 Broken power-law model 

The next continuum model we tried is based on a broken 
power law (bknpl) instead of a power law. The fit of the 
bknpl model gives already a reasonable fit with a reduced 
X 2 of 1.20, which is similar to the value for the fit with a 
pl+diskbb model. The break energy is at 2.9 keV, with a 
photon index of 1.7 before the break, and 1.6 after the break. 
The latter photon index is similar to the photon index of the 
power-law component in the pl+diskbb model. 

A Laor component added to this continuum model 



gives a reduced x 2 of 1.14, where the break energy and pho- 
ton indices are consistent with those in the fit without the 
line. We note that the bknpl+Laor model fits the spec- 
tra equally well as the PL+Laor+diskbb model does, with 
the parameters of the Laor component being consistent be- 
tween the two models. The normalization of the line is more 
than 3<r different from zero. 

To see whether a soft thermal component is required, 
we included a diskbb component to the model. The fit 
with a bknpl+Laor+diskbb model gives a reduced x 2 of 
1.12, which is a marginal improvement of the fit without 
the diskbb component. In contrast with the Laor com- 
ponent, the diskbb component is not significant, with a 
normalization that is < 3<r different from zero. The nor- 
malization of the line is 6a different from zero for the 
bknpl+Laor+diskbb model. The line in the best-fit model 
(i.e. bknpl+Laor) has an equivalent width of 67.5 eV. 

We note that the part of the power law before the break 
at ~3 keV possibly affects the contribution of the disk black- 
body, and that therefore the diskbb component in the bro- 
ken power-law model is not as significant compared with the 
continuum model based on a power law without a break. 



3.5 Comptonization model 

We then replaced the bknpl component with a model that 
describes the hard emission as Comptonization of soft pho- 
tons in a hot plasma (comptt; Titarchuk 1994). Using a 
disk geometry for the Comptonizing material, the comptt 
model gives a poor fit bf/v = 753/500). A spherical ge- 
ometry does not provide a good fit either. Adding a Laor 
component gives a reduced x 2 01 1-34, with the normaliza- 
tion of the line being more than 7a different from zero. 

The COMPTT+Laor+diskbb model gives a reduced x 2 
of 1.15, with an inner disk temperature of ~0.4 keV and a 
diskbb normalization that is > 3a different from zero. The 
line in the best-fit model has an equivalent width of 59.9 eV, 
with a line normalization that is > 4a different from zero. We 
note that the seed photon temperature of the COMPTT com- 
ponent becomes more than 10 times smaller than in the case 
when no diskbb is added to the continuum. The 0.6-10 keV 
flux of this three-component model is 3.9xl0 -10 erg cm" 2 
s _1 , whereas the flux of the diskbb component is 8.9x 10~ 12 
erg cm -2 s _1 . 

3.6 Reflection models: PEXRAV and PEXRIV 

The presence of a broadened Fe-emission line between 6 and 
7 keV suggests that the X-ray emission may partly be due to 
reflection of hard X-rays by a relatively cool disk. We there- 
fore explored fits to the data that include this reflection. We 
first used the pexrav and pexriv models, which describe 
the Compton reflection by neutral and ionized material, re- 
spectively (Magdziarz & Zdziarski 1995). 

The incident spectrum is assumed to be a power law 
with a high-energy cut off. Initially, we allowed this cut- 
off energy to be free to vary, but it always converged to 
values much higher than 200 keV, outside the energy range 
covered by our instruments, and we therefore used no cut off. 
Further, we fixed all abundances to solar and the reflection 
factor to 1. Compared to pexrav, the pexriv model has 
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Figure 2. XMM-Newton PN (0.3-10 keV), MOS2 (0.3-9.4 keV), RXTE PCA (3-20 keV), and HEXTE (20-100 keV) spectra of Swift 
J1753. In the lower panels the best fit residuals of a pl+diskbb model (left) and a PL+LAOR+DISKBB model (right) are shown. The Laor 
line profile is fitted to the broad emission line that is apparent in the spectra at around 7 keV. 



two extra parameters: The ionization parameter £, which 
we left free to vary, and the disk temperature, kT<n s k, which 
we fixed to the default value of 3xl0 4 K. 

The fit with the pexrav model gives a reduced \ 2 of 
1.38, and with the pexriv model the reduced \ 2 = 1-30. The 
pexrav/pexriv models include only the effect of bound-free 
transitions in the reflected spectrum. To account for the iron 
emission line at around 7 keV (caused by bound-bound tran- 
sitions), we added a separate LAOR component to the reflec- 
tion models. The spectral fits improved significantly, with 
a reduced \ 2 °f 1-19 an d 1-18 for the pexriv+Laor and 
pexrav+Laor models, respectively. The normalization of 
the line is 9<r different from zero in the fit with pexrav, and 
8<t different from zero in the fit with pexriv. The equivalent 
width of the line is 181 eV and 187 eV in the pexrav+Laor 
and pexriv+Laor models, respectively. 

Adding a diskbb component to the reflection plus 
line model gives a reduced \ 2 of 1-14 for both mod- 
els, with a best-fit disk temperature kT[ n — 0.3 keV and 
0.4 keV for the neutral and ionized models, respectively. 
The diskbb component is statistically not required in the 
pexrav+Laor+DISKBB model, with a diskbb normaliza- 
tion that is < 3<t different from zero. In contrast, for the 
pexriv +Laor-|-diskbb model the diskbb normalization is 
7cr different from zero. The flux of the diskbb component 
is 5.1x 10~ 12 erg cm -2 s" 1 and 5.5x 10" 12 erg cm -2 s" 1 
(0.6-10 keV) for the neutral and ionized reflection models, 
respectively. 



3.7 Ionized reflection model: REFLION 

Since pexrav and pexriv do not include the contribution of 
bound-bound transitions to the reflected spectrum, we also 
fit the data with a model that describes the reflection by an 
ionized, optically thick, illuminated atmosphere of constant 
density (reflion; Ross & Fabian 2005). The atmosphere 
is illuminated by a power-law spectrum with an exponential 
cut off at high energy (the cut-off energy is fixed at 300 keV, 
and hence outside our energy range). This model includes 
ionization states and transitions for the most important ions 
at energies above 1 eV, like O III-VIII, Fe VI-XXVI, and 
several others. We fixed the iron abundance to solar. Re- 



flion only provides the reflected emission, not the direct 
component, and we therefore used the reflion+pl model 
to fit the spectra, where we coupled the photon index of the 
illuminating power-law component in reflion to the photon 
index of the PL. 

The reflion+pl model gives a reduced \ 2 of 1.14. 
In this case we do not need to include a component to 
fit the line, since the line is part of the reflected emis- 
sion. The unabsorbed flux of this two-component model is 
3.9x 10~ 10 erg cm" 2 s" 1 (0.6-10 keV), where the flux of the 
power-law component is 3.0xl0 -10 erg cm -2 s _1 (0.6-10 
keV). The reflection ratio in the 0.6-10 keV energy range is 
0.3. Statistically, a diskbb component is not required to fit 
the data; adding this component gives a reduced y 2 °f 1-13 
with a best-fit temperature kTi n = 0.4 keV and a diskbb 
normalization which is less than 3a different from zero. The 
95% confidence limits for the diskbb component implies a 
disk flux < 6.1xl0~ 12 erg cm -2 s _1 (0.6-10 keV). 

To gain information about the disk radius, we convolved 
reflion with the Kerrconv component (Brenneman & 
Reynolds 2005). The Kerrconv kernel is applied to the 
reflected component to account for relativistic effects in the 
disk. We fixed the two emissivity indices to 3 (i.e. similar to 
the index that we found in the line profile) and the break ra- 
dius to a large value (this is equivalent to a single emissivity 
law for the whole disk). We further fixed the outer radius of 
the disk to its maximum value. We found that the spin of the 
black hole is not constrained by these data, and since we are 
testing for the possibility of a truncated disk at a distance 
that is much larger than the radius of the ISCO, we fixed 
the spin to its maximum value of 0.998 and left the inner 
radius of the disk free to vary. This allows the disk radius 
to float from ~1.23 R g to its maximum value. We checked 
that fixing the spin to zero did not change the results sig- 
nificantly. In summary, we allowed only two parameters of 
Kerrconv to be free: The inner radius of the disk, and the 
inclination of the disk with respect to the line of sight. The 
(reflion+pl) *Kerrconv gives a reduced \ 2 of 1.13, with 
an inner disk radius of 256 R g , an inclination angle consis- 
tent with 0° , and the rest of the parameters consistent with 
the model parameters without the Kerrconv kernel. 
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Figure 3. In the various continuum models fit to the XMM-Newton/RXTE spectra of Swift J1753, a broad emission line is apparent 
between 6 and 7 keV. Here, we show a zoom in of the PN residuals in the fit with the PL+LAOR+DISKBB model. To show the broad 
emission line, we have set the normalization of the line to zero, and applied an additional rebinning to the plot. 



3.8 The broad iron line at around 7 keV 

In previous subsections we included a Laor component to 
the continuum model to fit a broad line at around ~7 keV. 
The line is significantly required by the fit, as indicated by 
the fact that the normalization of the line is different from 
zero at a confidence level of 3-14a, depending on the model 
we used to fit the continuum. As shown in the right panel 
of Fig. [3] the line is extending from ~5 keV up to ~9 keV, 
with an equivalent width of ~60-187 eV, depending on the 
continuum model used to fit the spectrum. 

The best-fit parameters of the line depend on the model 
that we used to represent the underlying continuum emis- 
sion. For instance, for the different best-fit models in Tab.[T] 
the inner disk radius varies from 5.5 R g to 14.6 R g and the 
inclination appears to be constrained to high values, i > 85°. 
Such high inclination seems to be incompatible with the fact 
that no eclipses have been observed in the X-ray light curve 
of Swift J1753. A possible explanation for this high inclina- 
tion is that the broad line that we observe is a combination 
of lines from more than one ionization stage of iron. We 
tested this by adding three Laor profiles to the pl+diskbb 
model as described in i]3.2l with the energy of each line fixed 
at 6.4 keV, 6.7 keV and 6.97 keV, respectively, to account 
for possible emission from neutral and lowly ionized, He- like, 
and H-like iron, respectively. We further coupled all the re- 
maining parameters of these three lines, with only the line 
normalizations left free to vary. This multiple line model 
gives a reduced % 2 of 1.14, with the 95% confidence limits 
of the inclination angle 83° and 86° , and the line normaliza- 
tions of the 6.4 keV and 6.7 keV lines being consistent with 
zero. A more realistic approach would be to link the relative 
normalizations of the iron lines such that their ratios would 
be determined by the ionization balance in the disk. We note 
that the reflion component, as described in £13.71 calculates 
self-consistently the relative intensities of lines from different 
ionization stages of all elements in the reflected spectrum by 
fitting the ionization parameter of the disk. Recall that the 
fits of the reflion component convolved with Kerrconv 
in ^3.71 yield a large inner disk radius (~250 R g ) and a low 
inclination angle (i < 8° at 95% confidence level). 



The exact values of the line parameters may also de- 
pend on the approximations done to calculate the relativistic 
profile of the line. To test this, we tried three other compo- 
nents, besides the Laor profile, to represent the line profile: 
diskline (Fabian et al. 1989), Kerrdisk (Brenneman & 
Reynolds 2005), and Kyrline (Dovciak, Karas & Yaqoob 
2004). We only applied these components to the pl+diskbb 
model, such that we either fit pl+diskbb+diskline, 
pl+diskbb+Kerrdisk or pl+diskbb+Kyrline. For all 3 
models we left the inner radius of the disk, the inclination 
and the line energy free to vary, although the line energy 
was constrained to range between 6.4 and 7 keV. We fur- 
ther assumed a disk with a single emissivity index which 
we fixed to 3. For the models with Kerrdisk and Kyrline 
we fixed the spin of the black hole to 0.998. The results 
are unaffected if we fix the spin of the black hole to zero. 
As determined from these relativistic line profiles, the in- 
clination of the disk is 89.9° (with a 95% confidence lower 
limit of 69.7°), 90° (with a 95% confidence lower limit of 
89°), and 77.6°, (with the 95% confidence limits of 64.6° 
and 86.1°), for the diskline, Kerrdisk, and Kyrline pro- 
files, respectively. For comparison, for the Laor profile we 
found a 95% confidence lower limit for the inclination of 
~84° (see Tab.[TJ|. 

The value of the inner disk radius inferred from the line 
profile also depends somewhat on the model that we used 
to describe the line. For instance, the pl+diskbb+Laor 
model gives R[ n — 14.6 R g , with a 95% confidence upper 
limit of 19.7 R g (see Tab.QJ. For the pl+diskbb+diskline 
model we find Ri n = 11.7 R g , with a 95% confidence range 
of 6.2-15.5 R g , whereas the pl+diskbb+Kerrdisk and the 
pl+diskbb+Kyrline model give R ln = 1.9 R g , with a 95% 
confidence range of 1.7-3.8 R g , and Ri n — 15.5 R g , with a 
95% confidence range of 11.8-20.1 R g , respectively. 

We note that we used here the same data as those used 
by M06, and our analysis is similar to theirs; while here we 
report on a broad iron emission line with an equivalent width 
of 90 eV. Using the same continuum model, M06 reported 
an upper limit of 60 eV for the Fe line. We repeated the 
analysis of Swift J 1753 where we applied the same spectral 
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resolution used by M06 (a minimum of 10 counts per bin) 
and we simultaneously fit the PN, PCA and HEXTE spectra 
with a diskbb+pl model affected by absorption. We found 
best-fit parameters and reduced \ 2 values consistent with 
those in M06. We then added a Laor component to the 
model with the same parameter settings as described in H3.ll 
The diskbb+pl+Laor model gives x 2 /v = 2003/1951 with 
the line parameters consistent with those that we found. The 
line has a normalization that is more than 6a different from 
zero, and an equivalent width of 73.5 eV which is, within the 
confidence limits, consistent with the 90 eV measurement 
that we report. 



4 DISCUSSION 

We analyzed simultaneous observations of Swift J 1753 in 
the low/hard state taken with XMM-Newton and RXTE in 
March 2006. At variance with the results of Miller et al. 
(2006a), we can successfully fit the X-ray spectrum of Swift 
J1753 with continuum models that do not require emission 
from a prominent disk extending down to the innermost sta- 
ble circular orbit. Models including emission from a disk 
component, as proposed by M06, are also consistent with 
the data, although for those models where the disk compo- 
nent is required we found that the flux of the disk is a factor 
of 2 to 3 less than in the analysis of M06. Regardless of 
the continuum model, we found a broad iron line at around 
7 keV with an equivalent width of 60-187 eV, depending 
on the model. If the line is from a single ionization stage of 
iron and the broadening is due to relativistic effects near the 
black hole, then the line profile implies an inner disk radius 
between ~2 R g and ~16 R g , depending on the relativistic 
line profile (i.e. Laor, diskline, Kerrdisk or Kyrline), or the 
model used to fit the continuum. The inner disk radius can 
be larger if the width of the line results from the combination 
of lines from different ionization stages of iron, with each of 
these lines being broadened by relativistic effects. In fact, a 
model of reflection from an ionized disk that calculates the 
strength of all relevant lines self-consistently, convolved with 
a kernel that accounts for the relativistic effects close to the 
black hole yields an inner disk radius of 256 R g , with a 95% 
confidence lower limit of 246 R g . 

To get a better insight in the accretion disk geometry 
there are several factors that must be considered. We know 
that black-hole binaries show two basic states, the low/hard 
and the high/soft states (see e.g. Tananbaum et al. 1972), 
with transitions in between these states (see e.g., McClin- 
tock & Remillard 2006; Homan & Belloni 2005). The X-ray 
emission in the HSS is dominated by a thermal component 
which is usually identified as emission from the accretion 
disk. In this state, the inner radius of the disk, R[ n , ap- 
pears to remain constant despite large changes of the disk 
flux (e.g., Tanaka & Lewin 1995; Mendez, Belloni & van der 
Klis 1998). Tanaka (1992) interpreted this lack of change in 
_Ri n as a signature of the innermost stable orbit around a 
black hole. The black hole spectra in the LHS is dominated 
by emission from a power-law like component; in those cases 
where the interstellar absorption towards the source is low, a 
soft component with a low temperature (<0.5 keV) appears 
to be present in the spectrum (e.g., XTE J1118+480; Mc- 
Clintock et al. 2001). This component has been interpreted 



as due to a cool accretion disk truncated at a large radius 
(e.g., Chaty et al. 2003). This picture is strengthened by the 
fact that during transitions from the HSS to the LHS the ra- 
dius of the disk increases with a decreasing disk temperature 
(Gierliriski et al. 2008), suggesting that the disk recedes and 
cools down during the transition. Other sources show a simi- 
lar behavior (e.g., Kalemci et al. 2004), but since the inferred 
temperatures of the disk are below ~0.3-0.5 keV, interstel- 
lar absorption and inadequate coverage of the low-energy 
part of the spectrum make some of those results less secure. 
Changes of characteristic frequencies in the power density 
spectra of BHCs during state transitions (e.g., Mendez & 
van der Klis 1997; Homan & Belloni 2005) also seem to sup- 
port the idea of a receding disk, although in this case there 
is no clear-cut explanation about the origin of the variability 
(see the discussion in M06), and hence the evidence is less 
compelling. 

The picture of an accretion disk truncated at large radii 
in the LHS is also expected in some models. For instance, 
in the advection-dominated accretion-flow (ADAF) model in 
which the accretion flow is divided into two zones: A geomet- 
rically thin, optically thick accretion disk with a large inner 
radius, and a vertically extended inner region, the ADAF 
region, which is hot and optically thin and radiates less ef- 
ficiently (Narayan, McClintock & Yi 1996; Esin et al. 1997; 
see Narayan & McClintock 2008 for a review). This model 
has been successfully applied not only to sources in the LHS, 
but to quiescent sources as well (Narayan et al. 1996) . Some 
recent work (Meyer et al. 2007; Liu et al. 2007) suggests, 
however, that a part of the ADAF material may recondense 
back into a disk, but the contribution of this condensed disk 
to the total X-ray luminosity would be small (Taam et al. 
2008). 

In the last five years, XMM-Newton and Chandra re- 
vealed relatively broad iron emission lines of BHCs with pro- 
files that are consistent with being broadened by relativistic 
effects near the black hole (see Miller 2007 for a review). 
These broad iron lines are generally detected in the HSS 
(e.g., GRO J1655-40; Baluciriska-Church & Church 2000; 
XTE J1650-500; Miller et al. 2002b; GX 339-4; Miller et 
al. 2004; Dunn et al. 2008), in which the disk is thought 
to extend down to the ISCO (as described before), and the 
relativistic effects are therefore strong. Broad iron emission 
lines, consistent with a disk extending close to the ISCO, 
were also detected in GX 339—4 in the LHS (Miller et al. 
2006b), and in V4641 Sgr when the source was in an inter- 
mediate state (Miller et al. 2002a). We found a broad iron 
line in the XMM-Newton/ 'RXTE spectrum of Swift J1753, 
making it only the second case of a BHC in the LHS state 
showing such a broad line. A relativistic broadened iron line 
in this black hole state would be in contradiction with the 
standard scenario as described above, in which the disk is 
truncated at very large distances from the black hole. 

This contradiction appears to be resolved by recent 
studies of the BHCs Swift J1753, GX 339-4 and XTE 
J1817-330 (M06; Miller et al. 2006b; Rykoff et al. 2007) 
that reveal the possible existence of soft emission below ~2 
keV in these systems. This emission has been interpreted by 
these authors in terms of a prominent accretion disk in the 
LHS, with temperatures of 0.2-0.3 keV, and inner disk radii 
of ~2-6 R g . The idea of a cool disk extending all the way 
down to the ISCO in the LHS has recently been challenged 
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by Gierliriski et al. (2008) for the case of XTE J1817-330 
and by D'Angelo et al. (2008) for the cases of Swift J1753 
and GX 339—4. These authors pointed out that a truncated 
disk would be affected by irradiation from the same Comp- 
tonized photons that produce the high-energy emission, and 
would therefore appear to be hotter and a have smaller ra- 
dius than in the non-irradiated case. This effect was not 
taken into account in the fits of M06, Miller et al. (2006b), 
or Rykoff et al. (2007). Gierlihski et al. (2008) fit the Swift 
spectra of XTE J1817— 330 with a simplified model that ac- 
counts for irradiation, and they find that if this effect is 
not considered, the inner disk radius is underestimated by 
a factor of 2-3. D'Angelo et al. (2008) analyzed the case of 
an ADAF flow inside a truncated accretion disk and, from 
qualitative fits to the X-ray spectra of Swift J1753 and GX 
339—4, they conclude that reprocessing of hard photons by 
a truncated disk would mimic a soft disk extending down to 
very small radii. 

Using the same data as analyzed by M06, we found 
that there are realistic models that fit the X-ray spectrum 
of Swift J1753 without the need of a thermal disk-like com- 
ponent extending down to very small radii. We point out 
that for our fits we do not consider the effect of possible 
irradiation of the disk by high-energy Comptonized photons 
(Gierlihski et al. 2008; D'Angelo et al. 2008). We can calcu- 
late the magnitude of this effect by noting that in the case of 
the reflion+pl model only ~30% of the power-law flux is 
reflected off the accretion disk, and therefore the remaining 
~70% of this flux must be thermalized and re-emitted by 
the disk. From our results we deduced that the disk flux due 
to this irradiation would have to be ~2 xl0~ 10 erg cm~ 2 
s _1 (0.6-10 keV). However, in the fit with the reflion+pl 
model we did not detect emission from the accretion disk, 
with an upper limit that is a factor ~25 less than what is 
expected. (We note that the discrepancy is even larger given 
that one should also include the intrinsic -gravitational- disk 
flux; see Gierlihski et al. 2008.) This a challenge to the in- 
terpretation of the Fe emission line as due to reflection. 

In the case of GX 339-4, Tomsick et al. (2008) found 
that the disk component is significantly required to fit the 
LHS spectra at luminosities of 2.3% and 0.8% of the Edding- 
ton luminosity, a factor of ~2 and ~7 below the luminosities 
as probed by Miller et al. (2006b) for the same source. Tom- 
sick et al. (2008) found inner disk radii of ~4 R g however, 
following the results of Gierlihski et al. (2008) and D'Angelo 
et al. (2008), their values would probably increase by a fac- 
tor of ~2-3 if irradiation of the disk by the Comptonized 
emission was considered. Tomsick et al. (2008) also detected 
broad features due to iron Ka in the LHS of GX 339—4; their 
results indicate that if the width of the line is produced by 
relativistic effects in the disk, then the line must originate 
from within 10 R g of the black hole. 

In this work we used relativistic line profiles to fit the 
broad emission line at around 7 keV, although one should 
bear in mind that there are alternative models that ex- 
plain the width and the profile of the iron Ka line with- 
out the need of relativistic effects. For instance, Laurent & 
Titarchuk (2007) proposed a model that explain the proper- 
ties of the iron line in terms of down-scattering of hard pho- 
tons in a Comptonizing outflow with optical depth greater 
than 1. Done & Gierlihski (2006) applied a similar model to 
fit BeppoSAX data of the galactic BHC XTE J1650-500. 



They found that the line at around 7 keV in the spec- 
trum of a bright LHS observation can either be fit with a 
model including extreme relativistic effects, or by resonance 
iron K-line absorption from an outflowing disk wind and 
an emission line that is compatible with a truncated disc. 
On the other hand, Rozahska & Madej (2008) calculated 
atmosphere models for an accretion disk around a super- 
massive black hole irradiated by a hard X-ray power law, and 
they found that the observed spectrum contains a Compton 
shoulder that can contribute to the asymmetry and equiv- 
alent width of some observed Fe Ka lines in active galactic 
nuclei. One should also take into account that, as we have 
shown here, the best-fit profile of the line depends on the 
underlying continuum assumed to fit the data, which adds 
extra uncertainties to the inferred parameters of the disk. 



5 CONCLUSIONS 

We have shown that the X-ray spectrum of Swift J1753 can 
be fit with a continuum model that does not require a disk- 
like component, and inferred from the line profile the disk 
does not extend down to the ISCO, but it is consistent with 
a disk truncated at a few to a few hundred gravitational 
radii. However, as demonstrated by M06, the data of Swift 
J1753 also allow for fits with models in which the disk is 
truncated at radii close to the ISCO, although irradiation of 
the disk by the Comptonized emission may affect this picture 
(Gierlihski et al. 2008; D'Angelo et al. 2008). We found that 
in those cases that a disk-like component was required, the 
contribution of the disk to the total flux is a factor of 2 to 
3 times less than the flux contribution that M06 found for 
the disk component. 

We have detected a broad iron line in the spectrum of 
Swift J 1753. The inner radius of the disk deduced from fits to 
the line, using a Laor profile, seems to suggest a disk extend- 
ing down to 5.5-12.5 R g , for the best-fit models. Fits with 
different relativistic line profiles suggest somewhat larger in- 
ner disk radii, up to ~16 R g (diskline and Kyline), or lower 
Rin, down to ~2 R g (Kerrdisk), whereas reflection models 
smeared by relativistic effects suggest even much larger radii, 
Rin ~250 R g . The XMM-Newton/RXTE data as analyzed 
here, do not provide a definitive answer to the question of 
the accretion disk is truncated at large radii or not, since the 
answer strongly depends on the continuum and line model 
that are used. One should however bear in mind that the in- 
terpretation of the line profile in terms of relativistic effects 
close to the black hole is not the only possible one, and other 
alternatives do not require a disk extending down close to 
the ISCO to explain the shape of the line. It remains to be 
seen whether this kind of models can fit the data of Swift 
J1753 and other sources in the LHS. 
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Table 1. Results of various continuum models. In the first part of the table we show the results of a power-law (pl) based model. A line component (LAOR.), a disk blackbody (diskbb) or 
both were added (these components arc also added to the other continuum models). The PL parameters are the photon index, V, and the normalization JVp (in units of photons/keV/cm 2 /s 
at 1 kcV). The Laor parameters are the line energy, E^, which we constrain to range between 6.4 and 7 keV, the emissivity index, index, inner disk radius, R in (with R g =GM/c 2 ), 
inclination, i, and normalization of the line, TVl (in units of photons/cm 2 /s). The DISKBB parameters are the temperature at inner disk radius, fcTj n , and disk normalization, TVrj, which 
is defined as [.R m (km)/d(10 kpc)] 2 cosi, where d is the distance to the source. Next, we show the results of a broken power-law (bknpl) based model, with Ti, the photon index for 
energies below the break energy, break point, E^, photon index after the break, T2, and normalization, Ny (same units as PL normalization). Next, we give the results of a continuum 
model based on a Comptonization (comptt) component. The COMPTT parameters are: soft seed photon temperature, fcTo, plasma temperature, kT, optical depth of the plasma, t, and 
normalization, 7V C . Next, we show the results of two reflection models PEXRAV and pexriv, with parameters: photon index, T, cosine of the inclination angle, cosi, and normalization, 7Vp 
(same units as PL normalization). An extra parameters of the pexriv component is the disk ionization parameter, f. Finally, we give the results of another reflection model: reflion+pl. 
To gain information about the disk radius, we convolve this model with a Kerrconv component. The reflion components are the ionization parameter, and normalization of the 
reflected spectrum, TVr. The photon index of the REFLION component is coupled to the PL photon index. The free Kerrconv parameters are the inclination and inner radius of the disk. 
We fit each of the above models including an absorption model (phabs) with TVjj the equivalent hydrogen column, and a constant component to normalize between the instruments. 
For all continuum models, the unabsorbed flux is 3.9xl0 — 10 erg cm -2 s — 1 (0.6—10 keV). In the models with a DISKBB component, the flux of this component ranges between 0.3-1.6 X 
IO" 11 erg cm~ 2 s~ x (0.6-10 kcV). 



Model 



r 



(io- 2 ) 



(keV) 



. Laor _ 
-Rin 

(Rg) 



(deg) (10- 4 ) 



fcT in 
(keV) 



AT, 



D 



PHABS 

(10 21 cm" 2 ) 



PL+DISKBB 
PL+LAOR 
PL+LAOR +DISKBB 



1.60±0.01 5.57±0.05 
1.66±0.01 6.05±0.02 
1.62±0.01 B.71±0.04 



R ,„+0.16 
6 ' 40 -0.0 



1 36+ ' 31 
liSD -0.13 



90.0 
86.2 



0.0 
-1.8 
+0.1 
2.3 



5.8±0.7 
1.8±0.3 



0.38±0.04 36.7±10.4 
0.36±0.02 35.2^27 



1.61±0.02 
1.60±0.01 
1.63±0.02 



603.7/500 
646.1/497 
564.2/495 



Ti 



(keV) 



N v 
(IO" 2 ) 



El 

(keV) 



. Laor _ 
-Rin 
(Rg) 



(dog) (IO" 4 ) 



fcT in 
(kcV) 



. DISKBB . 



PHABS 

N H 

(10 21 cm" 2 ) 



BKNPL 

bknpl+Laor 
bknpl+Laor+diskbb 



1.69±0.01 2.9±0.2 
1.70±0.01 2.7±0.2 
1.71±0.04 2.9±0.2 



1.60±0.01 
1.62±0.01 
1.62±0.01 



6.17±0.04 
6.18±0.05 
6.36±0.34 



°-' a -0.12 
6 89+ ' 11 



o 7 +6.3 
' -0.9 
c+6.1 
-1.2 



3.97 



17.4±7.9 



+°-\ 1.8+O.i 



86.2 

86.2^3 



c+0.7 
J -0.4 



0.18±0.02 2458.8tt9og'.l 



1.67±0.02 

1.67±0.02 

2 09+ 16 
^■ ua _0.07 



601.8/500 
564.0/495 
552.9/493 



fcT 
(keV) 



kT 
(kcV) 



JV C 

(io- 2 ) 



(keV) 



index 



. Laor _ 
Rin 
(Rg) 



(deg) (IO" 4 ) 



fcT in 
(keV) 



PHABS 

N H 

(10 21 cm" 2 ) 



comptt+Laor 
comptt+Laor+diskbb 



0.13±0.01 

01+ 05 

u u -0.01 



11.9±2.4 



3.2±0.3 

^-0.03 



1.7±0.2 

4 1+ 01 

* -0.8 



6.79±0.20 
fi -,.+0.25 
b75 -0.35 



, „ + 2.5 
3- 5 -0.9 
, ,. + 6.5 
d - 5 -0.7 



12.4^ 



86.2iJ;i 
86.2+.°;* 



2.7±0.6 
1.6±0.6 



0.38±0.04 35.1+ 2 7 '^ 



1.55±0.03 
1.60±0.04 



662.8/495 
566.4/493 



. PEXRAV/PEXRIV . 



(erg cm/s) 



ATp 
(IO" 2 ) 



El 
(keV) 



. Laor _ 
Rin 
(Rg) 



(deg) (10~ 4 ) 



kT ia 
(kcV) 



PHABS 

n h 

(10 21 cm" 2 ) 



x 2 /" 



PEXRAV 


1.67±0.01 


0.16±0.04 




6.06+0.03 


















1.61+0.02 


693.7/501 


pexrav+Laor 


1.68±0.01 


0.13±0.05 




6.11+0.04 


6 40+ 12 


3 4+ 07 


5.5+0.8 


86.9±2.7 


3.4+0.6 






1.64+0.02 


590.0/496 


pexrav+Laor+diskbb 


1.65±0.01 


01+ 05 
u,u -0.0 




5.90+0.07 


6 40+ ' 12 


3.2+0.5 


5 6+ 1 - 5 


84.8 


+2.4 
-10.8 


2 5+ ' 6 


0.31+0.05 




1.68+0.06 


560.9/494 


PEXRIV 


1.68±0.01 


0.13±0.03 


134t 7 f 


6.09+0.03 


















1.64+0.02 


651.6/500 


pexriv+Laor 


1.69±0.01 


0.12±0.03 


0.0+0.01 


6.12+0.03 


6 40+ 1 


3.4+0.6 


5.5+0.7 


84.1 


,+ 1.4 
-2.3 


3.5+0.7 






1.65+0.02 


586.4/495 


pexriv+Laor+diskbb 


1.64±0.01 


o-oii°o:°o 5 


317+ 118 
^'-205 


5.85+0.07 


6 43+ - 57 
u '*°-0.03 


io-ol° 6 :°2 


5.6+3.2 


78.0 


+ 8.3 
-18.2 


1 0+ 09 


0.35+0.03 


si.oj;^ 1 


1.67+0.04 


562.2/493 



(io- 2 ) 



_ REFLION 

(erg cm/s) 



(10- 8 ) 



(deg) 



Rin 

(Rg) 



fcT in 
(keV) 



No 



PHABS 

N H 

(10 21 cm" 2 ) 



REFLIOK+PL 
REFLIOK+PL4 



1.54+0.01 4.10+0.15 
1.59+0.01 5.17+°lg 
1.55+0.01 4.18+0.26 



1.54+0.01 
1.59+0.01 
1.55+0.01 



505 3 



, + 528 
650 



7.6+0.8 
8.2+0.2 
7.9+0.3 



17 7 +18.8 
11 ■ '—11.0 



1.75+0.01 
1.66+0.01 
1.75+0.01 



568.4/500 
560.3/498 
562.1/498 



